Enhancement of thermal transport in the degenerate periodic Anderson model 
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The low-temperature transport coefficients of the degenerate periodic SU{N) Anderson model are 
calculated in the limit of infinite correlation between /electrons, within the framework of dynamical 
mean-field theory. We establish the Fermi liquid (FL) laws in the clean limit, taking into account the 
quasiparticle damping. The latter yields a reduced value of the Lorenz number in the Wiedemann- 
Franz law. Our results indicate that the renormalization of the thermal conductivity and of the 
Seebeck coefficient can lead to a substantial enhancement of the electronic thermoelectric figure-of- 
merit at low temperature. 

Using the FL laws we discuss the low-temperature anomalies that show up in the electrical resis- 
tance of the intermetallic compounds with Cerium and Ytterbium ions, when studied as a function 
of pressure. Our calculations explain the sharp maximum of the coefficient of the T -term of the 
electrical resistance and the rapid variation of residual resistance found in a number of Ce and Yb 
intermetallics at some critical pressure. 
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I. INTRODUCTION 



The low-temperature charge and thermal transport of 
heavy fermions and valence fluctuators with Ce, Eu, Yb 
and U ions display interesting and complex behavior, 
like a striking correlation^'^ between the low-temperature 
Seebeck coefficient a{T) and the specific heat coefficient 
7 — Cv /T. In many of these systems the dimensionless 
ratio q = |e| liniT^o c«/7^ is nearly the same, although 
the absolute values of 7 and a/T vary by orders of mag- 



nitude. In metallic systems, the data show; 



3.4.5.6.7 



small 



deviations from this universal behavior, due to the vari- 
ation in carrier concentration, and in bad metals, q can 
become quite largei*^. The Kadowaki- Woods (KW) ra- 
tio^ which is defined as p{T) / {'^T)'^ , where p{T) is the 
electrical resistivity, exhibits similar universal features, if 
one takes into account the effective low-temperature de- 
generacy of the /multiplet, as defined by the multiplicity 
of the crystal field (CF) ground state^iiiS, the carrier con- 
centration and the unit cell volumeii. 

The near constancy of the KW and q ratios brings to 
the fore the validity of the Wiedemann- Franz (WF) law, 
np/T = Co, and a possibility of enhancing the electronic 
thermoelectric figure-of-merit in strongly correlated ma- 
terials, ZT = a^T/np, where k is the electronic contri- 
bution to the thermal conductivity and Cq = 7r^fc^/3e^ 
the Fermi liquid (FL) Lorenz number. When the WF 
law holds, metals must have a thermopower larger than 
155 /iV/K to achieve ZT > 1; to date no metal has been 
found with so large a thermopower. In the temperature 
window where the effective Lorenz number (C — np/T) 
is reduced, one can achieve ZT > 1 with substantially 
lower thcrmopowers, which might make it possible to find 
strongly correlated metals that can be used for cooling 
applications at low temperature. 



The above-mentioned universality of the KW and q ra- 
tios, typical of a Fermi liquid (FL) state, attracted con- 
siderable theoretical attention. Kontani^ explained the 
KW ratio using the orbitally degenerate periodic An- 
derson model that he solved by the quasiparticle (QP) 
approximation of Yosida and YamadaF^ii^. Neglecting 
both the vertex corrections and the momentum depen- 
dence of the self-energy, he derived a KW ratio that 
depends on the degeneracy of the / states and brings 
the experimental data closer to the universal (theoret- 
ical) curved. Miyake and Kohnoi^ calculated the q ra- 
tio for the same QP dispersion as in Ref. [13, using an 
effective #-fold degenerate free-electron model in which 
the on-site correlation Uff is accounted for by the renor- 
malized hybridization between the c and / states. Re- 
stricting the average number of / electrons to nj < 1, 
they treated the QPs as free fermions and assumed that 
the repeated impurity scattering gives rise to an energy 
dependent relaxation rate. Neglecting the QP damping 
due to Coulomb repulsion, Miyake and Kohnc^** calcu- 
lated the low-temperature thermopower as a logarithmic 
derivative of the frequency-dependent conductivityi^ and 
show that the q ratio is quasi-universal number. 

The effect of electron-electron scattering on the trans- 
port coefficients has recently been studied by Grenze- 
bach et alJ^, using the dynamical mean-field theorjsii 
(DMFT) of the periodic spin-1/2 Anderson model. The 
auxiliary impurity problem generated by the DMFT was 
solved by the numerical renormalization group (NRG) 
method, which discretizes the energy spectrum and de- 
fines the temperature as the difference between the two 
lowest energy states. This provides accurate results for 
the static properties at arbitrary temperature but can- 
not provide the thermal transport in the FL regime. The 
DMFT-I-NRG calculations indicate an enhancement of 
the figure of merit when the temperature is reduced be- 



low its value at the resistivity maximum, and show a 
breakdown of the WF law due to electron correlations. 
Recently, we calculated the thermopower and the g-ratio 
of the periodic Anderson model taking into account the 
CF splitting^'* and enforcing ny < 1 at each lattice site 
(a realistic description of the CF states requires a lo- 
cal constraint). Solving the auxiliary impurity problem 
in the non-crossing approximation (NCA) we obtained a 
semi-quantitative description of the experimental data on 
heavy fermions and valence fluctuators in the incoherent 
regime. However, a detailed analysis of the DMT-I-NRG 
and the NCA results shows that nether method can de- 
scribe the thermal transport much below the Kondo tem- 
perature Tk nor establish the FL laws. 

In this paper, we discuss the coherent thermal trans- 
port of the periodic Anderson model with SU{N) sym- 
metry in the limit of an infinitely large Coulomb repulsion 
between / electrons. Such an effective A^-fold degener- 
ate model applies to intermetallic compounds with Ce, 
Eu, Yb and U ions in the FL regime, where the excited 
CF states can be neglected^^. For a given compound, 
the number of effective channels depends on the pressure 
and doping. We assume that the model has been solved 
in thermal equilibrium, so that the FL scale Tq ~ I/7 is 
known. The total particle number per unit cell n(^) is 
also n known; /i denotes the chemical potential. Alter- 
natively, we can assume that the values of 7 and n are 
taken from experiment. Using these equilibrium quanti- 
ties, and enforcing the constraint n^ < 1 at each lattice 
site, we construct an analytic solution for the stationary 
heat and charge transport in the FL regime. 

Unlike in previous work, we do not calculate the trans- 
port properties in the QP representation, because the 
operator algebra in the projected Hilbert space is not 
fermionic and the representation of the exact charge and 
heat current density operators is cumbersome. In the 
original fermionic representation, Mahan^S has demon- 
strated that the correlation functions between the above 
current density operators can be expressed in terms of 
transport integrals which differ only by powers of the 
excitation energy in their integrand, so that the same 
techniques can be used for their evaluation as in weakly 
correlated systems. 

The transport coefficients are obtained by the DMFT 
and expressed in terms of the average conduction electron 
velocity, the renormalized density of conduction states, 
and the frequency and temperature dependent relaxation 
rate which explicitly takes into account the QP damping. 
Since we are considering the FL regime, we use the QP 
approximation of Yosida and Yamada-- to relate all these 
quantities to the FL scale Tq and to the unrenormalized 
density of c states, A/J(w), evaluated at the shifted chem- 
ical potential fj,L = fJ- + A/i. The shift Afi is determined 
by the Luttinger theorem. Using the Sommerfeld expan- 
sion, we obtain the universal FL laws which show that 
the KW ratio depends not only on the multiplicity of the 
/ state and the average FS velocity, but also on the car- 
rier concentration and the unit cell volume, as observed 



experimentally—. As regards the g-ratio, we find that 
changes in the carrier concentration, induced by pressure 
or chemical pressure, lead to deviations from universal- 
ity. We also find large deviations from the WF law due 
to the lowering of the effective Lorenz number, which can 
lead to a substantial enhancement of the thermoelectric 
figure-of- merit ZT > 1. The change in the effective de- 
generacy of the model induced by pressure explains the 
pronounced maximum of the coefficient of the T^ term of 
the electrical resistance^ii^i^i^i and the rapid variation 
of residual resistanc o^^'^^ , found in a number of Ce and 
Yb intermetallics at some critical pressure. 

The usefulness of the approximate analytic solution of 
the DMFT equations becomes apparent when we real- 
ize that a numerical calculation of the above mentioned 
universal ratios encounters serious difficulty, and neither 
the NRG, nor exact diagonalization, nor quantum Monte 
Carlo approaches provide accurate transport coefficients 
in the FL regime (especially when the coherence temper- 
ature is low) . Combining analytical and numerical results 
enables a reliable estimate of transport coefficients at ar- 
bitrary temperatures, which is needed if the model is to 
be compared with experimental data. 

The rest of this paper is organized as follows: Sec. |TT] 
describes the DMFT calculations of transport coefficients 
in the low-temperature limit, in Sec. IIIII we use the QP 
approximation to calculate the renormalized density of 
states and transport relaxation time in the FL regime. 
Sec. IIVI establishes the FL laws, in Sec.|V]we use our re- 
sults to discuss the experimental data, and Sec. |Vl] pro- 
vides the summary and conclusions. 



II. FORMALISM FOR THE TRANSPORT 
COEFFICIENTS 

The 5'C/(7V)-symmetric periodic Anderson model is 
written in the standard form^^. 



H — Hb, 



and 



n 



I77ip 



■Hr. 



-^iJV. 



(1) 



Here, Tihand describes the conduction (c) band with un- 
perturbed dispersion Ck, assumed to be a even function 
of k. Because we have an SU{N) symmetric model, there 
are N distinct flavors of conduction electrons, which are 
labeled by index cr, so that 



N 



T-iband - 2J zJ ^kCko 



-k<T- 



(2) 



CT=1 k 



The operators c]^^ (cka) create (annihilate) a band elec- 
tron with momentum k and flavor a. The non-interacting 
density of conduction states JV^{uj) is calculated for each 
flavor and any Ck as A/'°(w) = l/(A/'jV) X^k ^(^ ^ ^k)' 
where A/i is the number of lattice sites and V the volume 
of the primitive unit cell. The characteristic bandwidth 
of the unperturbed c-DOS is W. The number of con- 
duction electrons is Nric (per unit cell), i.e., there are 



Tic conduction electrons of each flavor. We assume an 
infinite Coulomb repulsion, Uff —^ cxd, and describe the 
localized 4/ states by Himp, 



N 



w™p-EE^/^/-/-^' 



(3) 



where E'l labels the site energy of flavor a and V projects 
onto the subspace with zero or one total /-electrons per 
site. The number of /electrons is restricted to Uf < 1, 
i.e., there are rif/N localized electrons per flavor. The 
hybridization between same flavor c and / electrons is 
described by TCmix 



N 



Hr, 






where the conduction-electron operators are now written 
in real space. The total electron number operator isAf = 

'EaJ2ti4aCia + fLfia) and the chemical potential fi is 
adjusted to keep the total particle number n = Nuc + Uf 
constant (as a function of temperature or pressure). For a 
degenerate paramagnetic state, all flavors are equivalent 
and the label a can be dropped. 

The intermetallic compounds with 4/ ions are de- 
scribed by the iV-fold degenerate model with on-site hy- 
bridization, Vij — VSij^ which simplifies the calculations. 
For large iV, an exact solution is provided by the mean- 
field or the slave boson solutionis but for small N the 
correlations give rise to the Kondo effect which cannot 
be described by mean field theories. The physically rel- 
evant value of N depends on the effective degeneracy of 
the 4/ state. A single electron (or hole) in the lowest 
spin-orbit state of the 4/ shell of Ce (Yb) is 2 J -I- 1-fold 
degenerate, where J = 5/2 for Ce and J = 7/2 for Yb. 
This degeneracy is further reduced by CF splitting and 
for small hybridization, and temperatures at which the 
excited CF states are unoccupied, the value of N is de- 
termined by the degeneracy of the lowest CF state. If 
the hybridization is such that the CF splitting does not 
occur, or the excited states are thermally occupied, the 
effective degeneracy is defined by the lowest spin-orbit 
multiplet, which is a sextet for Ce and an octet for Yb. 
In all the physically relevant cases, 2 < iV < 8, the sys- 
tem exhibits the Kondo effect. 

At high temperatures, the qualitative solution can be 
obtained by perturbative scalingi^ which shows that the 
properties depend on an exponentially small energy scale 
Tk — W^exp(— 1/iVg), where g is the dimensionless cou- 
pling constant g = V'^Afc{fi)/Ef. At temperature Tk the 
perturbation theory breaks down. The low-temperature 
thermodynamic quantities of the SU{N) model are given 
very accurately by the DMFT-I-NRG and we assume 
that we know the numerical value of Fermi liquid scale 
To — n'^kg/SjV, where the linear coefficient of the low- 
temperature specific heat, 7, is measured per unit vol- 
ume. Alternatively, Tq can be estimated from the slave 



boson solution or taken from experiment. Unlike the sit- 
uation for the static thermodynamic quantities, no cur- 
rently available numerical method provides a reliable so- 
lution for the transport coefficients at low temperatures. 
In this paper, we evaluate the heat and charge transport 
of the SU{N) symmetric Anderson model by the Fermi 
liquid theory ^^ which holds for T <Tq. 

For the static and uniform transport we need the q — > 
component of the charge and heat current density oper- 
ators which are obtained by commuting the Hamiltonian 
with the appropriate polarization operators^^. The total 
current density obtained in such a way satisfies 



eN 
lim T 



E'^kc[^(^)'=kc7(T'): 



(5) 



where Vk — Vck/^ is the band velocity of the non- 
interacting (and unhybridized) c-states. For constant 
hybridization the energy current density is given by a 
similar expression which shows that the model satisfies 
the Jonson-Mahan theorem^^. This allows us to express 
the charge conductivity by a{T) — e^NLu, the ther- 
mopower by a(T)\e\T = — L12/L11, and the electronic 
contribution to the thermal conductivity by k(T)T — 
N{L22 — -^12/^11)- I'^ each of these expressions we have 
introduced the (single-flavor) transport integrals: 



Lmn — I duj \ — — 



771+n — 2 



A(^,T). 



(6) 



where /(w) = 1/[1 -I- exp(/3a')] is the Fermi-Dirac distri- 
bution function, lo is measured with respect to the chem- 
ical potential /i, and A(ll},T) is a function calculated in 
Appendix |X] by the Kubo linear response theory. 

At low temperature {—df/duj) approaches a delta func- 
tion and the main contribution to the integrals in Eq. Q 
comes from the low-energy excitations within the Fermi 
window, |a;| < fc^T, and with the wave vectors in the 
vicinity of the Fermi surface. In the Fermi liquid state 
A{u!,T) can be calculated in a straightforward way (for 
details see Appendix [X)) which yields in the u!,T -^ 
limit 



A(c^,T) = i4AA,(w)T(w,T), 



(7) 



for a three-dimensional system. Here, Vp denotes the 
average of v^. over the renormalized Fermi surface of hy- 
bridized states, t{uj,T) is the transport relaxation time, 
given by the momentum-independent self-energy of the c 
electrons. 



T(a;,r) 



Iml],(a;+,r)^ 



(8) 



and Ncioj) is the renormalized DOS of c electrons (for a 
single flavor), 



-^^H-^E^^(k^' 



(9) 



The spectral function of the c electrons 



AciKuj) 



1_ 



-Im Gc(k,cj+) 



(10) 



is obtained by evaluating the retarded Green's function 
Gc(k;, z) just above the real axis, where z = w"*" = to + iS 
and (5-^0 from positive values. 

The Green's functions of the periodic Anderson model 
are defined by the time-ordered thermal averages in the 
usual way and are cast in the Dyson form using the 
equation of motion or the diagrammatic expansion. The 
Dyson equations for c and / Green's function rea d^^'^'^ 



Gc(k,z) 



z- Ef -Y.f{z) + ^i 



[z - Ek + Ai] [z-Ef- E/ (z) +fi]-V^' 



(11) 



and 



G/(k,z) 



Ek + M 



[z-ek + ^i][z-Ef-J:fiz)+fl]-V^ 



(12) 

where z denotes a variable in the complex energy plane 
and S/(z) the local self-energy of / electrons which de- 
scribes the renormalization of the (hybridized) / states 
due to local Coulomb interaction. The Dyson equation 
for the conduction electrons can also be written in the 
alternative form 



Ge(k,z) 



Z- Ek + M- Sc(z)' 



(13) 



where the local self-energy of the conduction electrons 
satisfies 



Sc(^) = 



V' 



z-Ef + ^i~T.f{z)' 



(14) 



This self-energy describes the renormalization of the 
unhybridized c states due to the scattering on the / 
states, i.e., Ec includes the hybridization and the corre- 
lation effects. The DMFT condition ensures S/(w+,T) 
to be local but the locality of Yjc{u)'^ ,T) also requires 
the hybridization to be momentum independent. The 
dc transport relaxation time is obtained by substituting 
T.c{uj+,T) into Eq. dH) which yields 



t{lo, T) ^{uj-Ef- ReS/ + nf + {5- ImEf) 
h ~ V^{S~Im E/) 



-, (15) 



where the T, w-dependence is due to S/(a;,T). To 
account for the residual impurity scattering, which is 
present in any sample, we assume that 6 in Eq. (jl5[) 
has a small but finite value. For temperatures such 
that |Im 'Sf{uj,T)\ ^ S this additional scattering is ne- 
glected and we discuss the temperature dependence of 
transport coefficients by setting (5 = 0. At T=0 the self 
energy vanishes and t{uj,T) becomes a constant. The 
pressure dependence of the residual resistance is dis- 
cussed using Eq. ([7]) with a constant relaxation time. 



which neglects the renormalization effects due to the re- 
peated scattering on the impurities. Where appropri- 
ate, the clean limit implies we are in the regime where 
d <^ \lm 'Sf{uj,T)\, whereas, we must examine the dirty 
limit, where 6 > \lm 'Ef{Lu,T)\, to describe residual resi- 
tance data. 



III. THE FERMI LIQUID APPROACH 

In the FL regime, the imaginary part of S/ is small 
and can be neglected when discussing the renormalized 
excitation spectrum in the limit T,lu — > 0. Expanding 
I]/(aj) through linear order in frequency we find in that 
limit 

uj - [Ef + Re I]/(a;) - ^] w (tj - u;f)ZJ^ + 0{uj^), (16) 



with Z 



f 



[I — d'S f / duj\uj=a] being the enhancement 
factor {0 < Zf < 1) and (i/ = [Ef + RcE/(0) ~ fi]Zf. 
The parameters Cbf and Zf appear in the the description 
of the model at low-temperatures. 

If we ignore the imaginary part of the self-energy, both 
Green's functions in Eqs. pTjl and P^ develop poles^ 
which define the QP excitations with wave vector k. Us- 
ing Eq. p^ we can write the secular equation for the QP 
excitations in terms of the renormalized quantities, as 



{uj-ek + fi)iuj-Ojf)-V^ =0, 



(17) 



with V 
roots to 



n^. 



the renormalized hybridization. The 
describe two QP branches 




(ek - M + ^f) ± V (^k - M - <^/)^ + 41^^ 



(18) 
separated by the hybridization gap 2V. These excita- 
tions are only defined for energies close to the chemical 
potential, where Im "Ef can be neglected. For large to the 
QP approximation breaks down and Eq. (J17p is unphys- 
ical. At the FS we have fl^ = or fl^ — 0, depend- 
ing on the value of Nric + Uf. For Ce and Eu systems 
we place the bare / level below the chemical potential, 
which puts fj, close to the top of the lower QP branch. 



From Q^ 



we see that ki^ is then close to the zone 



boundary, where e-^p—^ ::iW{W denotes the half- width 
of the unperturbed c band). For Yb compounds, the bare 
level is above /x, so that the lower QP branch is full and 
the chemical potential is close to the bottom of the upper 
branch. The equation Q^ = gives kp close to the zone 
center, such that ek^ — M — —W. Close to the FS, the 
QP dispersion is very weak and fl^ describes two heavy 
QP bands with a half-width equal to Vf/W. 

The above derivation treats the QP as a non- 
interacting Fermi gas with effective parameters ujf and 
Zf (or V). These parameters can be related to the lin- 
ear coefficient of the specific heat, which is assumed to 



be known. Using the expression for total QP density of 
states 



^''''(^^-VKrT.S^^-^ 



(19) 



and the fact that the QP excitations defined by Eq. pT]) 
are infinitely long-lived (Im S/(w) = at T = 0) we 
write 



TT^fcl 



J: 



6 



NJ^Q^iO) = 



1 



3V fcsTo' 



(20) 



where the A^-fold degeneracy of the system has been 
taken into account. The relationship between 7 and 
Af^^ is derived assuming that thermally excitated QPs 
increase the average energy per unit volume by approx- 
imately AE ~ {kBT)'^NN'QP{0). The FL parameter Tq 
introduced in Eq. pO| defines the low-temperature scal- 
ing behavior of all physical quantities. In many systems, 
it also sets the temperature at which the QP descrip- 
tion breaks down. This temperature might or might not 
coincide with the high-temperature Kondo scale of the 
system. By definition, the inverse of Tq is given by the 
specific heat coefficient or, equivalently, by the density 
of the QP excitations of the SU{N) model, which can 
be calculated very accurately in thermal equilibrium. As 
shown below, the low-temperature thermal transport is 
also characterized by Tq. 

The FL scale Tq can be related to the partial densities 
of / and c statesi^ by expanding Eqs. (fTTI) and P^ for 
small uj. This yields the spectral functions. 



and 



Ae(k,w) « al{uj)S{uj - n^) 



Af{k,u) « ai{u)6{uj~n^), 



(21) 



(22) 



where we used the fact that uj = fl^^ can only be satisfied 



with one of the roots, and introduced the coefficients 

-I -1 



and 



aliuj) 



Jf,.,\ — 



1 



V' 



{uj - Uff 



aii^) = Zf- 



V^ 



rak(w)- 



(23) 



(24) 



(25) 



When we substitute Eq. ([25)) into Eq. (|T9|) and use 
Eq. ([20]) . the summation of Ac(k, w) and A/(k, w) over 
the first Brillouin zone yields 



These coefficients satisfy the FL sum rulei^. 



a£(^) + Zj'aiiu) - 1. 



2 1 



V NkeTo 



--MM + ZJ^Mf{Q), 



(26) 



which shows that the FL scale is set by the product of 
the large enhancement factor 1/Zj and Aff{0). 

To relate To to the renormalized parameters Qj and 
Zf, we express, first, the renormalized c and /DOS in 
terms of the density of unhybridized band states. Using 
the identity d{Lu — r^j^ ) = S{ek ~ ek^)/|d^i(. /^^k| and 
(iJ7j/(iek = '^k('-^)L=a±' where Ck^ is the solution of 
Eq. P?)) for a given (small) w, we obtain from Eqs. (PT|) 
and (121) the result 



Mc{uj)=Ml \uj + [i 



V^ 



UJf 



and 



J^fiu 



ZfV' 

[uj - J)/]2 



Mc{uj) 



(27) 



(28) 



The expressions in Eqs. (|27p an d (1281) are similar in spirit 
and in form to those in Refs. [20 and [23, where the 
fully interacting DOS of the periodic Anderson model 
is equated to the DOS of a non-interacting U — 
model. The Green's functions of that model are defined 
by Eqs. ^ and ^ for S/(tj) = 0. Here, we find it 
more convenient to relate Mc{uj) to the DOS of an unhy- 
bridized conduction band which is obtained from Eq. (|13p 
for T,c{uj) — 0. At low temperatures, where only the 
hybridized bands have physical relevance, this auxiliary 
conduction band is just a convenient mathematical con- 
struct. But at high temperatures, where the model has to 
account for the scattering of conduction states on local- 
ized paramagnetic /states, the unhybridized conduction 
band is physically relevant. 

The auxiliary model provides the renormalized Fermi 
surface (FS) of the periodic Anderson model simply by 
the Luttinger theorem, i.e., from the fundamental Fermi 
liquid relation which states that the volume of the FS 
cannot be changed by interactions. Using the standard 
form of the c electron Green's function [in Eq. ([1^)1 ] and 
the fact that Ec(z) is momentum independent (this holds 
for momentum-independent hybridization in infinite di- 
mensions) we immediately learn that the FS of the aux- 
iliary band model (with n = Nric + Uf electrons per 
unit cell) coincides with the FS of the periodic Ander- 
son model with the same number of electrons per cell. 
The shape of the renormalized FS is obtained by solving 
Eq. (J17p a,t UJ = 0, which yields the implicit equation 



ekp = M + -^ 



(29) 



in terms of the non-interacting dispersion. We recall that 
/i is fixed by the condition 

n = iVric + Uf = NV I duj [7V;(w) + 7V/(w)]. (30) 
>.' — 00 

The auxiliary model has exactly the same FS, determined 
by the equation 



Eki? — fJ-L, 



(31) 



where /ii is obtained from the integral 



(32) 



i.e., ^L is the chemical potential of a conduction band 
with n = ric + rif electrons. The same chemical potential 
^ (^l) appears in Eqs. ^ and ^ [Eqs. ^ and (l32l) ]. 
because Luttinger's theorem ensures that the number of 
k-points enclosed by the FS coincides with the total num- 
ber of electrons in the system. Equations (|29p and (|3T|) 
give the shift 



AfJ,^ flL- /i 






(33) 



Substituting Eqs. ^7^ and ^^ into Eq. ^^ and using 
Eq. ([33)1 to eliminate V^/Cjf, we find the relationship 
between J)/ and Tq, 



ujf = Afi 



NM°UiL)VkBTo/2 ^ JV 
1 - NMO{fiL)VkBTo/2 - 2 



AfiAf°{^iLWkBn, 

(34) 

where the last expression neglects the exponentially small 
term in the denominator. The sign of uJf is set by A/i 
which is positive for Ce and Eu compounds and negative 
for Yb compounds. The enhancement factor is obtained 
from Eq. ^ as Zf ^ Afi Cof/V^. 

The ratio 2a)/ /Tq ~ A^VA/J(/iL) depends, for a given 
n{fi), on the renormalized chemical potential /i, the aux- 
iliary one fiL, and the unit cell volume V. While /i^ 
is easily obtained from n{fi)^ the value of fi is difficult 
to find without numerical calculations. Except, for very 
small 1/, large N, and for n/ ~ 1 and ric ~ 1/2, when 
each conduction band is close to half-filling. In that case, 
the renormalization does not change much the partial oc- 
cupancies of the / and c states, and we can approximate 
ric — V^^ dijjN^{Lo)^ which yields 



NV / dojM'^iuj) 



Nuc = Uf, 



(35) 



where we have implicitly assumed jjll > ^^^ which ap- 
plies to Ce and Eu compounds. A similar result can be 
obtained for Yb compounds after an electron-hole trans- 
formation. For large N the integral is small and since the 
integrand has a maximum around jj, (the auxiliary band 
is close to half- filling), it follows that A/i is small. Using 
J\f^{^L) — A/J(/i) we estimate 



V 



dujAf°{uj) ~ AfiVAf^ifiL) 



N 



(36) 



Thus, for large N and n/ ~ 1, we have a simple relation 
ujf ~ nfkBTo/2. For small N, and/or small n/, Eq. ([55]) 
still holds but we cannot claim that A/x is small. liAf^{uj) 
decreases rapidly for a; > /i, as it usually does, we can 
conclude A/iV7V°(/iL) < n//iV < A/iV7V°(/i) but can- 
not express A/iVA/J (/xl) in terms of n/ or relate ojj to 
To in a simple way, as in the large- A'^ limit. 



The auxiliary model with unhybridized conduction 
electrons has a simple physical interpretation for sys- 
tems with small hybridization and nj ~ 1, i.e., for heavy 
fermions with very low Kondo temperature. In such sys- 
tems the 4/ ions are in a well defined valence state, the 
low-energy dynamics is dominated by spin fluctuations, 
and Uf is temperature-independentji^. At low temper- 
atures, the renormalized FS is defined by the set of k 
vectors satisfying Ck = fi + A^. This FS is large, because 
it encloses (ric+nf/N) states of each flavor. At high tem- 
peratures, the Fermi surface is small, because it encloses 
only He conduction states (/ states are localized and do 
not contribute to the Fermi volume). Ignoring the fact 
that thermal fluctuations remove the discontinuity in the 
Fermi distribution function, we approximate the renor- 
malized FS by the FS of unhybridized conduction states. 
(This holds for T ^ W.) On the other hand, the condi- 
tion ek = /i defines a set of k vectors which are close to 
the FS of the auxiliary band with ric electrons. This FS is 
also small, because it encloses only ric states of each fla- 
vor. Since n/ and ric are temperature-independent, the 
shift A/i provides the difference between the FS of the 
high-temperature paramagnetic phase with Nric conduc- 
tion electrons (and n/ localized states), and the FS of the 
low-temperature FL phase with A^ bands containing each 
Tie + nf /N hybridized states. The above considerations 
show that the formation of the coherent QP bands is ac- 
companied by a 'jump' of the Fermi volume from a small 
to a large value, and that the large FS encloses Uf/N 
additional states. The auxiliary model approximates the 
conduction band of the periodic Anderson model and al- 
lows us to estimate this jump. 

So far, we neglected the imaginary part of self-energy, 
since we considered only the low-energy excitations at 
T = In order to calculate the transport properties at 
low but nonzero temperature, we have to include the QP 
damping which is given by the imaginary part of the self 
energy (in the clean limit). To estimate Im T,f{uj,T) at 
low temperatures, we use the diagrammatic analysis of 
Yosida and Yamadai^. In infinite dimensions, the ex- 
pansion is in terms of the local Green's functions, and 
produces the FL expression^iii^ 



-Iml]/(a;,r) c -[cu^ + {7rkBTfKN^l)[VMfiO)fr}, 

(37) 
where Tf is the irreducible 4-point scattering vertex for 
electrons with different flavors. Eq. ([57]) is a straight- 
forward generalization of the result produced by sec- 
ond order perturbation theory in which the bare inter- 
action Uff is replaced by the scattering vertex Tf. In 
the limit of large correlations, when the charge fluc- 
tuations are suppressed, the Ward identityi^"^ gives 
ZJ^ = (Ar-l)VA0^(O)r/,suchthat7= (7r2fc|/6)A^(A^- 

l)V[J\ff{0)]'^Tf, where we neglect the much smaller con- 
duction electron contribution to 7. Substituting into 



Eq. 1371) yields in the T, w ^ limit^ 



Im I]/(w,T) ~ — 27r 



(TrfcsT)^ 



(N ~ l)N^VJ^f{0){kBToy 



(38) 



Expanding Ec(w,r) in Eq. (|14p into a power series 
[with the real and imaginary part of S^(w,T) given by 
Eqs. mi) and (f38|) , respectively] we obtain the self en- 
ergy of c electrons in the FL form. That is, in the FL 
regime, the real part of Tic{lo,T) is a constant and the 
imaginary part is a quadratic function of T and lo. 
Inserting S]j(a;,T) in Eq. (fTS]) for t{uj,T) and using 

the leading 



term 



:w^/Arc(0) yields for r,t 



t{l,,T) 



where 



h{N^ 



l)N^VN^{f,L)k%T§ 



(uj - Hjff 



2'KCb'j 
r,{T){l-^f{l- 



7r2fc|r2 



{^keTf 



ro{T) 



h{N~l)N^VM^(,lXL)T^ 



27r3 



y2 ' 



(39) 



(40) 



The second line in Eq. (|39p emphasizes the fact that the 
limit w ^ is taken before T ^ and that E/(w+, T) is 
only known up to the w^ terms, so that the Sommerfeld 
expansion cannot be extended beyond the second order 
(up to that order both forms produce the same result). 
The first (second) bracket in Eq. (|39p is due to the real 
(imaginary) part of E/(w+,T). Unlike Ec(w,r), the FL 
laws produced by the Sommerfeld expansion of transport 
coefficients are affected not only by the w-dependence of 
Im E/(u;,r) but of Re S/(w,T) as wefi. 

In order to make an estimate of the slope of the renor- 
malized / DOS, which is needed for thermal transport, 
we invoke the DMFT condition. This condition also pro- 
vides a physical interpretation of the low-energy parame- 
ter Cjf. In DMFT, we compute the local Green's function 
from the local self-energy 



Gf{z) 



E 

k 



= V 



G/(k,z) 



(41) 



z- Ef+^i-T.f{z) - 



Z — i + fl 



map it onto the Green's function of an effective single 
impurity Anderson model with a hybridization function 

A(z) 



F(z) 



1 



-i?/+/i-A(z)-E/(z)' 



(42) 



and adjust A(z) to make Gf{z) and F{z) identical. The 
DMFT procedure works because the large-dimensional 
limit guarantees that the functional relationship between 
the local self-energy and the local Green's function for the 
lattice is identical to the functional relationship between 



the self-energy for the impurity and the impurity Green's 
function. In the limit where T — and w — + 0, we 
approximate A(tj) = iAq with the constant Aq < and 
write the DMFT condition at low frequencies as 



Affiuj) = --Im F(uj+) 



1 



A^ 



T^\Ao\ (u - cof)^ + Ay 



(43) 



where A/ = ^oZf- This Kondo-like form of Mf{uj) 
holds only for uj <^ ujf^ just like the quasiparticle dis- 
persion makes sense only for a; <C Tq. It cannot be ex- 
trapolated to higher frequencies, where the approxima- 
tion A(a;) — zAq does not hold. The approximate form 
given by expression Eq. (|43p has a maximum at ujf , where 
the exact / DOS has a gap; the exact DMFT spectral 
function has a maximum between jj, and a)/. The width 
of effective Kondo resonance is A^:. In the FL regime, 
UJ < ksT <C a)/, the transport coefficients depend on 
the characteristic energy scales defined by this effective 
Kondo resonance. 

The width A/, like ujf, can be related to the FL scale 
To of the periodic Anderson model. Equating the /-DOS 
and the effective single impurity DOS at w = yields 



7rA/.VAA,"(ML) = - 



(44) 



where x — ujf/Af. Solving for x produces the result 
27rA^VAA0(/ii) 



A/ = 



l±Vl-(27rA/iVAAO(/.tL) 



-.UJf 



(45) 



and we choose the negative sign. Since A/iVA/"" (^l ) was 
shown to be very small [see the discussion around Eqs. 
(f35l) and p6|) ] the root in the above expression can be 
expanded to produce the lowest order result. 



Af = 



"^f 



:^NkBTo. 



(46) 



The above expression, which follows from the DMFT con- 
dition, ensures that that the initial slope of the renormal 
ized / DOS is very small. 



dbJ 



^ ~ nAfj.VM'Jif^L) « 1. 
A/ 



(47) 



Finally, we remark that the FS average of the unrenor- 
malized velocity squared, v^., can be found from the in- 
tegral 



= f d'^kSik-kp) vl, 



(48) 



where the (S-function restricts the integral to the renor- 
malized FS. In infinite dimensions Vp is a constant for 
all fillings but in lower dimensions the change of the 
Fermi volume with pressure or temperature can mod- 
ify Vp and affect thermal transport. Numerical calcula- 
tions for the 3-dimensional periodic Anderson model with 
nearest-neighbor hopping on a simple cubic lattice give 
Vp — {tai/hyv^ where v'^ ~ 1.4 for n ~ 1/2 and v^ <$^ 1 
for n ~ 1. 
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IV. THE FERMI LIQUID LAWS 

The transport coefRcients of the periodic SU{N) An- 
derson model with infinite correlation between / elec- 
trons are obtained in the FL regime by making the low- 
est order Sommerfeld expansion of transport integrals, 
Eq. ([H]), and expressing the integrand A(w) in terms of 
the known expressions for v^^, Nc{oj), and t{u!,T). The 
algebra is straightforward (for details see Appendix IBJ) . 
and yields the transport coefficients as simple powers of 
reduced temperature T JTq. 

The specific resistance of N parallel channels obtained 
in such a way is 



P{T) = 



Ott^V 



he-^vlN{N - l)[NV^f^{^lL)Y 



tJ («' 



which holds for A^ > 2 and arbitrary A/'°(a;). 

The resistivity expression in Ea. p5)) neglects terms of 
the order of (T/Tq)'* and is only valid for T < Tq- Even in 
that temperature range, p{T) deviates from the universal 
(KW) form, because the pre-factor of the (T/Tq)'^ term 
has an explicit parameter dependence. The value of these 
parameters depends o the Fermi volume of the system, 
which can be changed in several ways. For example, addi- 
tional impurity scattering or lattice expansion (negative 
pressure) can localize the / states and exclude them from 
the Fermi volume. Another possibility is to increase the 
effective degeneracy of the / states (by pressure or ther- 
mal population of the excited states), which changes the 
number of resonant channels and shifts the Fermi surface 
away closer to the zone center. The 'jumps' in the Fermi 
volume changes fiL and Vp, which has a strong impact 
on the resistivity. This feature can be used to explain the 
resistivity anomalies that accompany the localization or 
delocalization of / electrons in heavy fermions, i.e., the 
breakdown or the formation of the QP bands. The rapid 
change of the coefficient of the T^-term in the resistivity 
following the pressure- or doping-induced 'jump' of the 
Fermi volume in Cerium compounds, is discussed in more 
detail below. On the other hand, if we tune the model 
parameters of the SU{N) model in such a way that the 
Fermi volume is preserved (by keeping N/n constant), 
the Luttinger theorem ensures that /ii and Vp do not 
change. In that case, the pre-factor of (T/Tq)^ is con- 
stant but p(T) changes due to variations in To which can 
be exponentially fast. 

For large A^, the approximation NVN^Xhl) — n^/A/i 
yields the expression 



P{T) 



Sl{A^x/nffV^ 



hN{N ■ 



-{iTf 



(50) 



which simplifies the discussion of heavy fermions in the 
Kondo limit, n/ ~ 1. The Luttinger theorem ensures 
n//A/i is constant even when the tuning of the exter- 
nal parameters gives rise to a charge transfer between 
the / and c states and thereby changes the renormalized 
chemical potential. As long as the tuning does not affect 



the degeneracy of the ground state, the pre-factor of the 
(7T)^-term of the resistivity is always the same. 

The Seebeck coefficient is obtained by using the Som- 
merfeld expansion for L12 and Ln which gives 



a{T) 



T 



3\NA^lV^foi^lL)To^ 



(51) 



The enhancement of L12 is solely due to the real part 
of the self energy. The imaginary part corrects Ln, and 
gives a factor (3/2) enhancement of a{T) with respect 
to the U = case. This factor does not arise in mean 
field theory or in slave boson approximations based on 
quadratic Hamiltonians with renormalized parametersi^. 
The above result corrects the expression used in our pre- 
vious paper^^, which was derived neglecting the energy 
dependence of the density of states. The term A^, which 
does not occur in the low-temperature resistivity expres- 
sion in Eq. (|49p . has an explicit parameter dependence, 
so that, strictly speaking, a{T) is not a universal func- 
tion of T/Tq. In bad metals and systems with a low- 
carrier concentration, fi^ is close to the band edge, where 
Af^{fiL) could be very small, so that a{T) could be very 
large. 

The Seebeck coefficient of heavy fermions with n/ ~ 1 
and large A^ assumes the simple form 



a(T) = T 



n/|e| To' 



(52) 



Since the doubly occupied / states are removed from the 
Hilbert space, the model is highly asymmetric, and the 
initial slope Iuht^o a(T)/T never vanishes. As a matter 
of fact, the closer the system is to half- filling, the larger 
is the slope, a/T oc I/Tq. This, however, does not nec- 
essarily imply a large thermopower, since the FL laws 
are only valid for T <C To and close to half-filling To is 
exponentially small. 

The FL result for the thermal conductivity in the clean 
limit reads 



where the usual Lorentz number, Cq 
been replaced by the effective one, 



Co{T) = Co 



32^2 



(53) 
7r^fc^/3e^, has 



(54) 



and Co — n kg/2e . This change is due to the imaginary 
part of the self energy and therefore not obtainable by 
mean-field or slave-boson calculations which neglect the 
QP damping. The T ^ limit yields the Wiedemann- 
Franz law, k{T) oc Ta{T), but the correction given by 
the square bracket leads to deviations even at low tem- 
peratures. Since the factor multiplying the T^ term is 
very large, we find a reduction of k{T) and substantial 
deviations from the WF law much below To. 



Before cfosing this sectfon, we summarize the proce- 
dure for calculating the transport coefhcients of the pe- 
riodic Anderson model in the FL regime. The model is 
specified by the hybridization V, the bare /-level position 
Ef, the degeneracy N, the center of the bare conduction 
band Eq, the unperturbed dispersion ek, and the restric- 
tion nf < 1. The unperturbed density of conduction 
states N^{iLj) is easily obtained from ek- For a given to- 
tal electron density n (per cell) we use the DMFT-I-NRG 
or some simpler scheme, like the slave bosons, to find the 
renormalized chemical potential, the number of / elec- 
trons, and the low-temperature coefficient of the specific 
heat which sets the FL scale To- These static quanti- 
ties can be calculated in thermal equilibrium to a very 
high accuracy. For constant hybridization, we can make 
separate DMFT-I-NRG runs for nc{fJ,) and nf{fi) and 
do not have to calculate the spectral function which is 
not given very accurately (for all cu) by the NRG. From 
n = Nile + nf, we obtain /i, /i/,, and A/i, which provide 
the renormalized FS and the average vj^. From Tq, /i^, 
A/x we obtain the Kondo scale (D/ and specify completely 
the low-energy behavior of Nc{io) and r(w). The FL laws 
follow at once. For heavy fermions with large N , Uf "^ 1 
and small V , the procedure simplifies considerably, since 
we can approximate AfiNVJ^^{nL) — nf, and obtain all 
the renormalized quantities using ujf ~ nfkBTo/2. 



KW ratio, can deviate from the universal valueii^. In bad 
metals and systems with a low-carrier concentration, fj.^ 
is close to the band edge, where M^{^l) could be very 
small, making a{T) large. For a given system, the q-ratio 
can be pressure dependent due to the transfer of / elec- 
trons into conduction band. The data on the pressure 
dependence are not available but the deviations from the 
universal value are indicated by recent chemical pressure 
data^. 

A further remarkable consequence of correlations is 
the enhancement of the low-temperature figure-of-merit 
due to the deviations from the WF law. Using Eq. ([55]) 
and neglecting phonons in the FL regime we express the 
figure-of-merit as a ratio ZT = a'^(T)/C{T). For con- 
stant Lorenz number, the maximum of ZT is defined by 
the thermopower but in correlated systems the tempera- 
ture dependence of the effective Lorenz number can lead 
to an additional enhancement. Even though our FL re- 
sult is valid only for T ^ Tq, it captures the essential 
features: an increase of a^ and a decrease of C{T) that 
ultimately give rise to ZT > 1. The enhancement of ZT 
is due to the renormalization of the thermopower and the 
Lorentz number i.e., ZT > 1 is not restricted to metallic 
systems with a{T) > 155 fiV/K. We expect large ZT for 
small To but to find the optimal situation one should tune 
the parameters and study the border of the FL regime 
by numerical methods. 



DISCUSSION OF EXPERIMENTAL 
RESULTS 

A. The universal ratios 



The FL laws in Eqs. (|49l - [53|) describe coherent charge 
and heat transport in stoichiometric compounds, in a way 
analogous to the phase-shift expressions for dilute Kondo 
alloys^. They explain the near-universal behavior of the 
KW ratioSiiS reported for many heavy xfermions and va- 
lence fluctuators. The ratio p{T)/{'jT)'^ obtained from 
Eq. (|^5)) exhibits an explicit dependence on the ground 
state degeneracy and the average velocity (squared) , and 
an implicit dependence on the Fermi volume, i.e., on 
the carrier concentration. The A'^-dependence as well as 
the dependence on carrier concentration of the KW ratio 
were recently emphasized by Kontani and his collabora- 
tors^' ^°, who obtained the power law n~^'^ for the latter 
by using the free electron approximation for the average 
velocity and the density of states in the expression for 
the resistivity. Our formulation leads to a similar A^- 
dependence [see also Eq. ([50)1 ] but includes also the de- 
pendence on the average velocity, the renormalized Fermi 
volume and the carrier concentration, which are implicit 
functions of the degeneracy and have to be taken into 
account when discussing the pressure or doping experi- 
ments on strongly correlated electron systems — . 

As regards the thermal transport, Eq. ([?T|l gives q — 
limT^o \ea\/-fT ~ 12/[NAfjM^{nL)], which has an ex- 
plicit parameter dependence, so that the g-ratio, like the 



B. Pressure dependence of the low-temperature 
resistivity 

The FL laws derived for the periodic Anderson model 
explain the changes observed in the transport coefficients 
of heavy fermions under applied pressure. As an il- 
lustration, we consider the pressure dependence of the 
coefficient of the T^ term in the resistivity, defined as 
A{p) = {p — pq)/T^, where po is the residual resistiv- 
ity. In the case of the two heavy fermion antiferromag- 
nets CePd2-a;Ge2-a?^ and CeRu2Ge2^, A{p) is small 
and nearly pressure independent for p < 4 GPa. Above 
4 GPa, A{p) increases rapidly and reaches a maximum 
value for pressure between 4 and 8 and GPa. At the 
critical pressure, Pc, the maximum of A{pc) is typically 
one order of magnitude higher than at ambient pressure. 
For p > Pc, the value of A(p) drops to the ambient pres- 
sure one. A somewhat different behavior is found in the 
heavy fermion superconductor CeCu2Si2 and in the an- 
tiferromagnet CeCu2Ge2, when analyzed in the normal 
state^i,;^!. The values of A{p) are large at initial pres- 
sure which is applied in order to restore the FL behavior; 
above that pressure A{p) decreases to a plateau and then 
drops by nearly two orders of magnitud o^^'^^ . 

A sharp maximum of A{jp) is observed for systems in 
which the ground state of the 4/ ion is characterized at 
ambient pressure by a CF doublet well separated from 
the excited CF states. The Neel temperature of such 
systems is much higher than the Kondo temperature and 
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large paramagnetic entropy of Ce ions is removed at Tjv 
by an AFM transition rather than by a Kondo crossover. 
The low-entropy state involves large unscreened moments 
which correspond to the 4/ ions frozen in the high- 
temperature (magnetic) configuration. Below Tjv, the 
conduction electrons are essentially free, except for some 
magnon and impurity scattering, such that the resistiv- 
ity at ambient pressure is small and weakly temperature- 
dependent, as shown by CePd2-2;Ge2-2; and CeRu2Ge2 
data. 

To explain the experimental data we assume that an 
increase of pressure leads at pc to the delocalization of 
the 4/ states^* and the formation of hybridized bands, 
which can be described by the SU{2) Anderson model. 
The two channels (sub-bands) in which the lowest 4/ 
doublet hybridizes with the conduction states of the ap- 
propriate symmetry, accommodate nearly one electron 
per site, ric -I- n//2 ~ 1, for each flavor. The renormal- 
ized Fermi volume is large and the FS is close to the 
edge of the Brillouin zone, where Vp is small. The cor- 
responding value of fiL is of the order of the bandwidth 
[see Eqs. ((3T|) and ((32|) ]. such that Af^{fiL) is also small. 
[The maximum of A/'°(ti') is assumed to be close to the 
center of the band.] Since Tq is small (Tq < Tjv), the 
coefficient of the T^ term in the resistivity, given by Eq. 
(|49)) for N = 2, is large. The enhancement of A{p) with 
respect to the values at ambient pressure (where the / 
states are localized) is due to the delocalization of the / 
states by applied pressure. In a system with delocalized 
/ states, the main effect of pressure is to increase the hy- 
bridization and To, provided the degeneracy of the lowest 
occupied CF level is preserved. We can understand the 
increase of To by recalling that the FL scale of the lattice 
is proportional to the low-energy scale of the auxiliary 
impurity model [see Eq. (|34p ] and that an increase of the 
Kondo scale with hybrdization (pressure) is a typical fea- 
ture of any Kondo system. An increases of Tq can also 
be inferred from the experimental dat a^^'^^ , which show 
that in CePd2-sGe2-s and CeRu2Ge2 the value of Tq 
scales with the high-temperature Kondo scale. 

A qualitative change occurs at the point where the 
hybridization becomes so large that the system cannot 
sustain the CF excitations and the degeneracy of the / 
level changes from two to six. To estimate the resistiv- 
ity of this high-pressure state we use the SU{6) model in 
which a single / electron is distributed over six equiva- 
lent hybridized channels. In that case, our FL solution 
shows that the FS is shifted away from the zone bound- 
ary, ^L is decreased and N^{^l) increased, with respect 
to the values obtained for A^ = 2. The average squared 
velocity v^p and the FL scale Tq are also increased for 
N = 6, which reduces A{p) to small values, in agreement 
with the experimental data^ii^^iS^i^. The drop of A{p) 
for p > Pc signifies the doublet-sextet crossover and is 
mainly due to the pre- factor in Eq. (|49)) . which has an 
explicit and implicit dependence on N . At the crossover, 
the Fermi volume changes and the KW ratio is strongly 
pressure dependent. Once the degeneracy of the / state 



is stabilized at a higher value, a further increase of pres- 
sure reduces A{p) by increasing Tg but does not change 
v'jp or ^L, which are fixed by the Luttinger theorem for 
{ric + nf/d) states per channel. Eventually, at very high 
pressure, the system is transformed into a valence fluc- 
tuator with an enormously enhanced FL scale, such that 
A{p) drops to a very small value. 

As regards CeCu2Si2^ and CeCu2Ge2^, we can ex- 
plain the data using the same reasoning as above, if we 
assume that the / states are delocalized at the initial 
pressure which restores the FL state. The value of A{p) 
at the initial pressure is large, because the temperature 
dependence of the resistivity is due to two degenerate 
sub-bands with very heavy fermions. The FS of these 
sub-bands is close to the zone boundary, where Vp, ^l, 
and N'cif^L) are very small. The hybridized bands which 
involve the excited CF states have a small Fermi vol- 
ume and their low-energy properties are free-electron like. 
They carry most of the current but have temperature- 
independent resistivity and do not affect A{p). The first 
decrease of A{p) occurs at pressure at which an increase 
of the hybridization removes the CF excitations, such 
that the effective degeneracy of the / state increases from 
2 to 4 or 6. For large N, the expression in Eq. ^^ gives 
small A. The final drop of A{p) is due to the crossover 
into the valence fluctuating regime and an exponential in- 
crease of the FL scale Tq. If we assume a proportionality 
between Tq and the Kondo scale Tk, which is indicated 
by the experimental data, our calculations would explain 
the A{p) versus Tk scaling reported in Ref. |23. 

The above scenario also offers an explanation for the 
peak observed in the residual resistivity of CeCu2Ge2^ 
and CeCu2Si2^ for pressures at which the degeneracy of 
the / level changes and A{p) drops from the maximum 
value. At such a pressure, the formerly non-resonant 
channels responsible for the low residual resistivity trans- 
form into resonant ones, which reduces Vp and ^^^(/xl). 
For a constant impurity scattering rate, this leads to 
a drastic increase of po- A further increase of pres- 
sure transforms the system into a valence fluctuator with 
nearly free conduction states and small pQ. An alter- 
native explanation for this peak has been suggested by 
Miyake and Maebashi^S, who attribute it to critical va- 
lence fluctuations. In our approach, these would occur at 
higher pressure, and we expect the resulting feature to be 
narrower than the one observed experimentally. A peak 
in po{p) is also observed in YbCu2Si2^. In this case, 
the mechanism that leads to it is the "mirror-image" of 
the preceding one^: starting from a valence fluctuating 
regime at ambient pressure, the system is driven into a 
state with a well defined valence in which all eight com- 
ponents of the J = 7/2 4/-multiplet hybridize with con- 
duction states, so that po is very large. A further increase 
of pressure stabilizes the magnetic 4/^^ configuration^l, 
inhibits the charge fluctuations, and gives rise to the CF 
excitations. This reduces the effective degeneracy of the 
/ hole by splitting-off the doubly-degenerate 'resonant' 
sub-band with heavy QPs from the "non-resonant" sub- 
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bands with nearly free electrons. For the reasons given 
above, pq is reduced yet larger than at ambient pressure. 



VI. SUMMARY AND CONCLUSIONS 

A detailed understanding of the thermoelectric prop- 
erties of the periodic Anderson model in the full temper- 
ature range and for arbitrary parameters is of a consid- 
erable interest, as it might facilitate the search for new 
thermoelectrics with a useful low-temperature figure-of- 
merit. The model with the SU{N) symmetry and an 
infinitely large Coulomb repulsion between / electrons 
captures the low-temperature features of the intermetal- 
lic compounds with Cc, Eu and Yb ions which exhibit 
often a large thermopower and might have a potential 
for applications. 

The thermoelectric response of the SU{N) model with 
n = Nnj + Tie particles is calculated in the FL regime 
assuming that the FL scale To is known. The value of 
To can be obtained from the numerical DMFT results 
for the linear coefficient of the specific heat or estimated 
from an analytic expression for the QP DOS given by 
the slave-boson approximation^^. The transport inte- 
grals are related to the current-current correlation func- 
tion and calculated by using the Luttinger theorem and 
the DMFT condition. The low-temperature transport co- 
efficients are obtained from the Sommerfeld expansion as 
power series in terms of the reduced temperature T/Tq. 
The coefficients in the expansion depend on the average 
conduction electron velocity, the unit cell volume, the 
effective degeneracy of the / state, the unrenormalized 
density of c states, A/'"(/iL), and the shift in the chemical 
potential A/z = /il — /i which also has a simple physical 
interpretation. The chemical potential fi corresponds to 
n interacting particles and for a flat low-energy disper- 
sion, /i is not much different from the chemical potential 
of Nuc conduction states decoupled from the / states. 
The FS of such a conduction band encloses Nuc points 
in the k-space and is considered to be 'small'. The FS of 
the Anderson model, determined by the Luttinger theo- 
rem, accommodates / electrons in addition to c electrons 
and is considered to be 'large'. For k-independent hy- 
bridization, this 'large' FS coincides with the FS of a free 
conduction band with n electrons and chemical potential 
jj,L. Obviously, the FS of a non- interacting conduction 
band with n electrons differs from the FS of Nric elec- 
trons, so that the shift A/i measures the 'jump' in the 
Fermi volume due to the hybridization. 

The parameter dependence of the coefficients multiply- 
ing various powers of reduced temperature corrects the 
simple scaling behavior and explains the deviations of the 
KW and the g-ratio from the universal constants. The 
FL law for the conductivity gives the KW ratio which 
depends on the multiplicity of the / state, the unit cell 
volume, the average FS velocity, and the carrier concen- 
tration, in agreement with the experimental data^l. The 
FL law for the thermopower gives the g-ratio which de- 



pends on the concentration of / electrons. From these 
results we conclude that, in some systems, pressure or 
chemical pressure can cause a substantial shift of the 'uni- 
versal ratios' from the common values. We also find that 
the quasiparticle damping leads to the break-down of the 
WF law due to the temperature dependence of the effec- 
tive Lorenz number. In the absence of the thermal cur- 
rent due to phonons, this would lead to a substantial en- 
hancement of the thermoelectric figure-of- merit ZT > 1. 
Assuming that pressure changes the hybridization and 
the effective degeneracy of the / state we explained the 
pronounced maximum of the coefficient of the T^ term of 
the electrical resistanc e^^'^^i^^'^"^ and the rapid variation 
of residual resistanceSli^ found in a number of Ce and 
Yb intermetallics at some critical pressure. 

Our results describe the main features of thermal 
transport that one finds in heavy fermions in the FL 
regime but do not indicate the temperature at which the 
FL description breaks down nor relate the FL scale to the 
Kondo scale, which characterizes the high-temperature 
behavior. A complete discussion should explain the low- 
and high-temperature behavior on equal footing, provide 
all the relevant energy scales of the system, and account 
for the change in the effective degeneracy of the low- 
energy states due to applied pressure, doping or temper- 
ature. The change of the effective degeneracy can have 
a dramatic effect on the transport coefficients and seems 
to be responsible for complicated thermoelectric response 
of intermetallic compounds with Ce and Yb ions. The 
realistic modeling should remove the SU{N) symmetry, 
take into account the excited CF states, consider the de- 
tails of the band structure, and/or include aditional in- 
teractions. Unfortunately, none of the presently available 
methods can solve such more realistic models, explain the 
crossovers between various physical states, and describe 
the behavior of the correlation functions in the full tem- 
perature and pressure range. 

A rough description of the high-temperature regime, 
which takes into account the CF splitting, can be ob- 
tained by assuming that the conduction electrons scat- 
ter incoherently on the 4/ ions. This reduces the lattice 
model with n electrons per unit cell to an effective Ander- 
son impurity model which can be solved very accurately 
by the NCA^^'^'^. The NCA solution leads to an exponen- 
tially small Kondo scale and explains the main features 
that one finds in the thermoelectric response of heavy 
fermions and valence fluctuators for T > Tk- In particu- 
lar, the NCA shows that the effective degeneracy of the 
model changes as the excited CF become thermally pop- 
ulated. Some justification for applying an effective single 
impurity model to the stoichiometric compounds is pro- 
vided by the fact that in the temperature range we are 
concerned with the resistivity of most heavy fermions is 
very large and the mean free path is not much longer than 
a few lattice spacings. Furthermore, the NCA solution 
for an effective two-fold degenerate impurity modeUSii^S 
agrees for T >Tk with recent DMFT-f NRG results for 
the two- fold degenerate periodic Anderson model^^. For 
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an A''- fold degenerate model the NCA gives a{T) and 
the power factor P{T) — a^(T)/p(T) with a pronounced 
maximum and C{T) with a shallow minimum around Tr-, 
which supports our previous conclusions regarding the 
enhancement of the thermoelectric figure-of- merit. The 
agreement between the perturbative solution of an effec- 
tive single impurity model and the experimental data in- 
dicates that the high-temperature state of heavy fermions 
and valence fluctuators can be represented by a nearly 
free conduction band which is weakly perturbed by lo- 
calized (paramagnetic) / states. Ignoring the thermal 
broadening of the Fermi distribution function we find 
that this FS is much smaller than the low-temperature 
one which must include the / electrons in order to sat- 
isfy the Luttinger theorem. The perturbative solution 
breaks down in the coherent regime, as it cannot describe 
the change in the Fermi volume. For a periodic model, 
the reduction of the paramagnetic entropy, the crossover 
from the high-temperature perturbative regime into the 
FL regime, and relationship between Tk and Tq can only 
be obtained by non-perturbative methods. 

The DMFT+NRG method describes, in principle, the 
crossover from the high-entropy state formed above Tk 
to the low-entropy FL ground state state. It provides 
the Kondo scale Tk which governs the high-temperature 
behavior and gives an accurate numerical estimate of the 
FL scale Tq and other thermodynamic quantities, like 
the number of particles or the chemical potential of the 
ground state. However, neither the DMFT-I-NRG nor 
the effective impurity approach provide a quantitative 
description of the transport properties of the periodic 
Anderson model for temperatures below Tr- or Tq. The 
usefulness of the FL approach is that it gives the low- 
temperature transport coefficients for a given Tq and al- 
lows us to obtain an overall description by extrapolating 
between the FL solution and the high-temperature one 
obtained by the NCA or the DMFT-^NRG methods. 

It would be interesting to perform an experimental 
study of the pressure-induced deviations of the g-ratio 
and the KW ratio from their universal values, using the 
universal behavior of the power factor or the effective 
Lorenz number as consistency checks. Since P{T) and 
C(T) require only transport measurements, they are well 
suited for pressure experiments. The above discussion 
makes clear that pressure experiments provide the most 
stringent test of the FL laws. 
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APPENDIX A: CORRELATION FUNCTION IN 
THE FERMI LIQUID REGIME 

To find A(a;, T) we use the Kubo formula for the static 
conductivity which provides Ln as the zero-frequency 
limit of the number current -number current correlation 
function^^ 



raf3 _ 



= Q/3/ 



ill = hm Re— Li, (i^). 



(Al) 



where 



■ali 



Liiim) 



VM 



rfTe-'^(T.ji(r)j;3(0)), (A2) 



i/i = 2'KkBTl is the Bosonic Matsubara frequency, the r- 
dependence of the operator is with respect to the Hamil- 
tonian in Eq. ([1]), the subscripts a and (3 denote the 
respective spatial indices of the number-current vector 
jc/e, and we must analytically continue L"f (w;) to the 
real axis L"\^{u) before taking the limit ly ^ 0. 

Substituting the definition of the particle current op- 
erator from Eq. ([5]) into Eq. (|A2[) for L"£ and taking the 
limit of infinite dimensions, in which the dressed correla- 
tion function is equal to the bare one^^i^, we obtain 



Lt?{m) 



^ ' dTe''^'^VvkaVk/3G,(k,r)Ge(k,-T). 



V^f^Jo 



k 



(A3) 

where Gc(k, t) is the imaginary time Green's function 
of the c electrons, which can be expressed as a Fourier 
series Gc(k, t) = fc^T^j^ exp(— zcL'„T)Gc(k, ia;„). Sub- 
stituting into Eq. (jA3[) and integrating over imaginary 
time provides the result 



-a(3 



Liiim) 



-n 






ksT} ^ y ^VkaVk;3 



Gc(k, zcj„)Gc(k, zcj„ + ii^i) 



(A4) 



which has to be analytically continued to the real axis. 
Since the Green's functions in Eq. (|A4p depend on k only 
through e(k) which is an even function of k, while Vka is 
odd, the summation over k vanishes for a ^ (3. We con- 
sider only isotropic systems, where L"f = Ln. To per- 
form the analytic continuation we follow closely Ref. l3a 
and obtain the result (evaluated explicitly for a three- 
dimensional system) 



Lii = 



X Re 



'h vj. 



VAf,3 



d. hm ^(-^-^(- + -^ 



1/^0 



Gciuj) - Gciuj + ly) GI{lj)-Gc{lo + v) 

V + Sc('^) — Sc(tJ -^ v) // + S*(cij) — I]c(w + ' 



(A5) 
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where we introduced the local Green's function 

G,M = ^Ge(k,c.), (A6) 

k 

and replaced the square of the a-component of velocity 
by its Fermi surface average Vp/d {d is the spatial dimen- 
sion, which we can take to be equal to 3 for real systems) . 
This step is justified, because the energy integrations in 
Eq. (|A5[) is restricted to a narrow interval around the 
Fermi energy, where the integrand is singular, so that 
the main contribution to the k summation comes from 
the k-points close to the Fermi surface; on the infinite- 
dimensional hypercubic lattice, the integral can be per- 
formed exactly and one finds that the average square 
velocity is equal to aft*'^ /Wfi^ when expressed in terms 
of the reduced nearest-neighbor hopping t* = 2vdt. 
We can now take the limit of i^ — > 0. Writing 



l^^O V + Hc\y^) — Sc(w + V) OLU 






and 



lim 



G*{uj)- Gc {uj + iy) _ ImGc (uj) 



produces our final result 



ill 



duj 



dfiio] 



duj 



A{lo,T), 



(A7) 



(A8) 



(A9) 



where A(w,r) is defined by 



A(^) 



/ ImGc(w) 



Re 





dG,{uj) 


duj 


1 


as,(w) 


dtj 



VA/'cBtt I fniEc(w) 

^(AlO) 

To estimate the relative importance of the two terms in 
Eq. (jAlOp we introduce the Hilbert transform of A/'c((jj), 

1 



rtciuj) = - 



TrVTV" 



ReGc(w). 



(All) 



In the FL regime, where Afc{oj) is (5- function like, the 
slope of Ti.c{uj) is very large; it is proportional to the c 
electron enhancement factor 



Re 



1 






On the other hand, Im dY,c{ui)/duj is small around w = 0, 
because ImEc(cij) is close to its maximum value. Using 



Im 



doj 



<z: 



we neglect [Itti 9S]c(a;)/9w]^ in the denominator of the 
second term for K{lo) in Eq. (jAlOp and approximate 



Re 



duj 



1 






z, 






-im 



duj 



duj 



This term is small with respect to A/'c(w)/ImI](w) which 
diverges in the limit T, w — > 0. Keeping only the singular 
term in Eq. (jAlOp we obtain the result of Eq. ([7]) used in 
the text. 



APPENDIX B: TRANSPORT COEFFICIENTS 

The transport coefficients are obtained by expanding 
the derivative of the Fermi function in a Sommerfeld ex- 
pansion, -df{io)/dLo = 5{u}) + {Tr^klTy6)[d^S{cj)/duj^], 
which gives the following result for the transport integrals 

Lrnn = [c.™+"-2A(c., T)] ^^^ (Bl) 

'"+"-2A(cj,r)] 



7r2fc|r2 



6 



52 



duj^ 



OJ 



For thermal transport we have to evaluate 

Ln - [A(c.,r)Uo + ^!^A"(a;,r)U.o, (B2) 



L12 — 
L22 = 



6 

^^^A'(c.,T)U.o, 



(B3) 

(B4) 



where 



A(c^,T) = -vlKiL0)T{LU,T) 



(B5) 



A'(c.,T) = -4[AA^Mr(c.,T)+^,Mr'(a;,r)] 



(B6) 



A"(c.,T) = -4[AA!'(c.)r(c.,r) 



+ 2AA^Hr'(c^,T)+AA,Hr"(a;,r)], (B7) 

and A' and A" denote the first and second derivatives 
with respect to lo. The first derivative of A/'c(w) is 



KH = ^^" ' "^ V ,(c.) + (" ' ^'^' mUco), (B8) 



and becomes in the lo -^ Q limit 



Kio) 



-2 






(B9) 



where the DMFT condition showed that the derivative 
of the /-electron DOS is small at a; = 0. For the second 
derivative, we obtain 



A/'c'(0) = 



2AA,(0) 



u} 



(BIO) 



using Mf{0)l ZfV^ = A/'c(0)/(i? and dropping the terms 
proportional to CjfN'AQ) and a>?A/'f'(0) which are expo- 
nentially small [see also Eq. ((47|) ]. 

The transport relaxation time given by the expression 
in Eq. ^ yields at tj = 



r'(0,T) = -2 



ro(T) 



r"(o,r) = 2To(r) 



-? 



{^keTY 



(Bll) 
(B12) 
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The transport integrals are now easy to find. The first 
one is obtained from Eq. (|B7[) . which gives 



Ln ~ -vlMMMT)^ {1+ 



+ 



7r2fc|T2 



A(0,T) 



6AA,(0) 



^} 



{nkBTY 



0{T^) 



(B13) 



The last term in the square bracket in the second line 
grows as 1/T^, while the first one is a (large) constant 
which we neglect at low enough temperature. Since 
A(0, T) ex l/T^, this approximation amounts to keeping 
the T^ terms in the electrical resistance and neglecting 
the T'' contribution. The correcting factor 2/3 in Eq. 
(jBlSp originates from the imaginary part of the self en- 
ergy and it is well known from the dilute alloy problem—; 
hence it arises only in the clean limit. Inserting the ex- 
pressions for the renornialized c-DOS [Eq. (|27|) ] and the 
relaxation time [Ea. (|40p ]. we obtain the dominant low- 
temperature contribution to in: 



ill 



n h(N - 



l)[Nv^f°,{^iL)f 

97r3V 



(B14) 



which yields for the static conductivity a{T) — Ne^Ln 



a{T) 



he^vlNjN ~l)[NVJ\f^{^jiL)Y 
97r3V 



f ) • (B15) 



and the resistivity p{T) = l/cr(r). 

The second transport integral is obtained from 
Eq. (|B6P which gives 



L 



12 



^A(0,T) 



MM t{o,t) 



-2Li 






T4ii 



njTo 



(B16) 



(B17) 



In the second equation we used ojf ~ ±nfkBTo/2, which 
holds for heavy fermions with nj^ ~ 1 and large N only. 

The thermal conductivity is defined by the expression, 



-(T) = N^ 



L22 



L12 



(B18) 



and follows from the previous results which give 
Lii = 2A(0,r)/3, Li2/iii ^ 2n^klTyCjf, L22 = 
7r2fc|T2A(0,T)/3, and ^22/^11 = n'^klT^/2. In the 
heavy fermion limit this yields the expression in Eq. ([53|) 
used in the main text. 
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